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Abstract 

O" 1 It is well known that the Schwarzschild solution describes the gravitational field 

$h I outside compact spherically symmetric mass distribution in General Relativity. In 

^ particular, it describes the gravitational field of a point particle. Nevertheless, what 

is the exact solution of Einstein's equations with <5-type source is not known. In 
the present paper, we give mathematical meaning of nonlinear Einstein's equations 
with 5-type source and find asymptotically flat static spherically symmetric solution 
in the class of generalized functions of combined system of equations of motion for 
metric and point massive particle. 
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1 Introduction 

In this article, we consider the classic problem: find the gravitational which is produced 
by a point massive particle. If particle is at rest, then the gravitational field is spherically 
symmetric. The spherically symmetric solution of the vacuum Einstein's equations is 
well known: it is the Schwarzschild solution pp. Therefore, it is often stated that the 
Schwarzschild solution describes gravitational field of a point particle. This statement is 
incorrect because there is no 5-type energy-momentum tensor corresponding to a particle 
on the right hand side of Einstein's equations. 

On the other hand, the solution of Einstein's equations outside point massive parti- 
cle must be isometric to Schwarzschild solution. Therefore the natural question arises: 
"Where is the ^-function ?". The answer turned out to be unexpected: 5-function cor- 
responds to infinite value of the Schwarzschild radial coordinate. Namely, in this paper, 
we prove that the Schwarzschild solution in isotropic coordinates is the solution of Ein- 
stein's equations in a topologically trivial space-time 1R 4 with 5-type source. The obtained 
space-time is geodesically complete. 

Severe mathematical difficulties arise during solution of this problem. A solution of 
Einstein's equations must be understood in a generalized sense after integration with test 
functions because the 5-function stands on the right hand side. But there is no multiplica- 
tion in the space of generalized functions and the question arises what is the mathematical 
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meaning of the left hand side of Einstein's equations which are nonlinear. Besides, locally 
nonsummable functions arise during solution of equations, and some functionals must be 
attributed to them. In other words, regularization is needed. In the preset paper, the 
solution of Einstein's equations is found in a generalized sense: the equations are satisfied 
after integration with test functions. We choose the usual space D(IR 3 ) of infinitely differ- 
entiate functions on R 3 with compact support as the space of test functions. Spherically 
symmetric solution is found in the conjugate space V(M?) in which we use the analytic 
regularization for spherically symmetric functionals. Though no multiplication is defined 
in V, the left hand side of Einstein's equations is well defined. 

The obtained solution turned out to be the well known Schwarzschild solution in 
isotropic coordinates. It appears after gluing together two exterior solutions of a black 
hole along the horizon. This solution is isometric to Einstein-Rosen bridge and is asymp- 
totically flat not only at large distances but near the particle itself where the curvature 
of space-time tends also to zero. At large distances, the gravitational field is attractive. 
Attraction decreases near the particle and changes to repulsion after crossing the horizon. 
This repulsion results in geodesic completeness near the particle. 

Attempts to interpret the Schwarzschild solution in terms of generalized functions were 
made earlier [2-6]. Our approach is essentially different. 



2 The main equation 

We consider three-dimensional Euclidean space IR 3 with Cartesian coordinates fi = 
1, 2, 3, in which the rotational group SO(3) acts in the usual way. Let the positive defi- 
nite Riemannian metric g IJLV {x) which defines Christoffel's symbols and the corresponding 
curvature be given on it. In our notations, the main geometric notions have the following 
form: 

r M / := -g pa (d^g ua + d v g^ - d^g^), (1) 

u p L up i^i/J- pp J- pp J- ue \ L vp L pt i V / 

Rppv P i (3) 

g^Rpu, (4) 



r^pvp 
Rpv 

R 



where T^ u p are Christoffel's symbols, R^p" is the curvature tensor, R^ u is the Ricci tensor, 
and R is the scalar curvature. 

Consider the covariant equation 

y/\g\R= -lGnM6(x), M = const, (5) 

where 

5(x) := ^x 1 )^ 2 )^ 3 ) 

is the three-dimensional 5-function. The volume element y/\g\ on the left hand side is 
needed for covariance because 5-function is the scalar density with respect to coordinate 
transformations. 

We shall see in what follows that equation (jSJ) is the main equation in the system 
of Einstein's and point particle equations. Moreover, this equation may have its own 
applications in differential geometry which are not discussed here. 
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We are seeking spherically symmetric solutions of equation (jSJ), and that is in agree- 
ment with the symmetry of the right hand side. Any spherically symmetric metric is Weyl 
flat, which means that there exists the coordinate system where the metric has the form 

where 8^ u := diag (+ + +) is the Euclidean metric and = <p(x). Then expressions for 
the Ricci tensor and scalar curvature follow from formulae fl2])-fl3J): 

R/iv = - d^dycf) + 5^ (A(f) + d<p 2 ) , 

where <9 M := d/dx^ is the partial derivative, A := 5^ u d^d y := d^d^ is the usual Laplacian 
in the Euclidean space K 3 , and we introduced notation dip 2 := S^d^d^. Using the 
expression for the scalar curvature, the Ricci tensor can be written in the form 

R P v = d^d u (j) - d^(j)d u (j) + ^<V<90 2 + -g^R. (7) 

The main equation fl5]) for the Weyl flat spherically symmetric metric takes the form 

4 A0 + 2 = -16irM5(x). (8) 

We rewrite this equation in another form which is suitable for the analysis. Introducing 
the new variable 

/ := e* (9) 

we obtain the following equation 

Af-^r = -4nM5(x). (10) 

Thus the main equation fl5]) for the Weyl Euclidean metric has the form ( FlOl) . It is 
nonlinear because of the second term on the left hand side. Solution of this equation 
should be understood in the generalized sense after integration with a test function (see, 
i.e. [7]) because there is the 5-function on the right hand side. If we look for a solution in 
the space of functionals (generalized functions) V then serious problems arise. Since the 
equation is nonlinear, the generalized functions must be multiplied, but multiplication in 
T>' is absent. Second, we shall see in what follows that some terms in equation ( 1101) are 
not locally summable, and hence the regularization of the appearing integrals is necessary. 
Though multiplication of functionals in T>' is not defined in general, for some particular 
functionals the left hand side has nevertheless definite meaning. 

We look for spherically symmetric solution of equation (Tl0|) in the spherical coordinate 
system r, 9, ip. In the space of spherically symmetric functionals r A , A G C, we use the 
analytic regularization [8j. Then the left hand side of the main equation (1101) is well 
defined for these functionals. 

Theorem 2.1. A spherically symmetric asymptotically flat solution of equation DD\) exists 
in the space of generalized functions T>' , for which the left hand side of equation / TiOj) is 
defined. It has the form 

M M 2 ( M\ 2 



r 4r 2 \ 2r / 
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Proof. Simple calculations lead to the following expression for the second term on the left 
hand side of equation ([ID]) 

df 2 f' 2 M 2 , s 

t — I — do) 

2/ 2f 2r 4 ' K 1 

where prime denotes differentiation with respect to radius. Therefore equation ( ITOj) for 
the generalized function (11 ip takes the form 

A / M M 2 \ M 2 e , . 

The function M/r is the unique fundamental solution of the Laplace equation decreas- 
ing at infinity (see., for example, [7]): 

/ M\ M 
A 1 + — = A— = -AtcMS(x). (13) 
\ r J r 

The function M 2 /2r 4 is not locally summable and requires regularization. We use 
the analytic regularization for the generalized function r A , A G C [8J. For re A > —3 the 
function r A is locally summable, and the functional 

(r A , ip) := / dxr X Lp 

is defined for all test functions ip G X>(1R 3 ). This functional is analytically continued on the 
whole complex plane A G C except simple poles located on the real line A = —3, —5, —7, . . . 
[8]. The poles does not matter in the considered case because we do not get into them. 
The equality 

Ar A = ^d r (r 2 d r r x ) = A(A + l)r A ~ 2 

is simply verified for re A > 2. It remains valid also after analytic continuation to the 
point A = —2. For A = —2, it looks as follows 

4 = ^- 

Thus, for analytic regularization, we attribute to the locally nonsummable function 1/r 4 
the functional 

The last functional is well defined. 

We see that the nonlinear term on the left hand side of the equation ( jTOi) is canceled 
after integration with a test function 

A M 2 M 2 

A ¥J " = °- (14) 

if the analytic regularization is used for the function M 2 /2r 4 . □ 
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Equation (jSJ) for a massive point particle was considered in [2-4] for Weyl Euclidean 
form of the metric. The authors used parameterization / := x 2 an d obtain equation 



X A X = -2nM5(x) (15) 
which is equivalent to equation ( TTUj) . They proposed the solution 

* (r) = 1+ 2^0) (16) 

differing from that corresponding to solution (fTTj) by the essential factor x(0). This factor 
is obtained from the regularized ^-function and depends on small parameter e, 

x(o)~y T , 

which is the "radius" of the regularized 5-function. In the limit e — > 0, solution f lT6|) 
becomes trivial, x = 1j an d satisfies equation y^i? = corresponding to a vacuum 
rather than point particle. So the solution proposed in [2-4] does not depend on mass 
M. A generalization of equation (jSJ) for a charged massive particle interacting with 
electromagnetic field was also considered. In this case, the factor x(0) becomes nontrivial, 
and the total mass of a particle is proportional to its charge. This effect is interpreted as 
the regularization of the self-energy of a point charge by gravitational interaction. 
To show the difference in the approaches, let us consider equation (|T5|) for 

*=! + £, (IT) 
which is different from ([TBI) . Substitution of this solution into equation (fT5|) yields 

1 + — J A ( 1 + — ^ = l 1 + — 1 (-2ttM)5(cc) = -27rM5(a;). (18) 
2r / \ 2r J \ 2r J 

Expression -S(x) is not defined. We define it as follows. Relation 

Ay = A(A- l)r x ~ 3 + r x A^ 

can be easily verified for re A > 3. We analytically continue it to the point A = — 1 where 
it can be rewritten as 

1.1 1 2 
-A- = A— - 



So, we define 



1 8(x) (p) — 1 ( 1 A 1 (p) — 1 ( A 1 2 p 

For analytically regularized \ the right hand side is zero, and we get relation 

-S(x) = 
r 
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which is valid in a generalized sense. In other words, the relation r x 5(x) = which is 
valid for A > is analytically continued to the point A = — 1. 

We see that solution (TTTj) can be obtained either by solving equation ( TTOj) or ([15!) . In 
both cases, we used the analytic regularization and get the same solution. Our solution 
fill I) does depend on mass M. It is nontrivial for neutral particle and differs essentially 
from solution (fTO]) considered in [2-4]. 

Thus we have found the generalized solution of equation (flO]) in the space of func- 
tional T>'(IR 3 ). For this solution the left hand side of the main equation is well defined. 
Unfortunately, we did not describe explicitly that subspace in D'(IR 3 ), for which the left 
hand side of the main equation is defined. This question is complicated and related to 
the uniqueness of the obtained solution. We leave it for future work. 



3 Relation to the Schwarzschild solution 

Consider four- dimensional topologically trivial space-time with coordinates t, p, 9, <p. We 
identify sections t = const with three-dimensional space which is considered in spherical 
coordinate system p, 9, <p. Static spherically symmetric solution of vacuum Einstein's 
equations without 5-type sources is well known. It is called the Schwarzschild solution pQ 
and is usually written in the form 

ds 2 = (l _ ^ _ d P^ M _ p 2 ( ^2 + sin 2^ 2) (lg) 



P 



The Schwarzschild solution is defined for 2M < p < oo (outside the horizon p s = 2M) and 
for < p < 2M (under the horizon). It is asymptotically flat and tends to the Lorentz 
metric when p — > oo. It is important that constant M entering the Schwarzschild solution 
is an integration constant and is not contained in the vacuum Einstein's equations at 
all. From mathematical point of view, constant M may take arbitrary values M G M. 
However, if we assume that the Schwarzschild solution describes the gravitational field 
outside a point particle, then comparison with the Newton gravitational law at large 
distances (see, i.e. 0) tells us that the integration constant M is the mass of a particle 
and therefore must be positive. 

The space part of the metric is not Weyl Euclidean, and hence the radial coordinate is 
denoted by the letter p. To compare it with the three-dimensional metric obtained in the 
previous section, we perform the coordinate transformation r i— > p(r) and change the total 
sign of the metric to make it positive definite. Then the space part of the Schwarzschild 
metric ffl9|) becomes 



J'2 



dl 2 = P 2M dr 2 + p 2 (d9 2 + sin 2 9d V 2 ). (20) 
1 



P 



To make the spatial part of the Schwarzschild metric Weyl Euclidean, the function p(r) 
must satisfy the equation 



P 12 t 
2M r 2 ' 
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For p > 2M, it reduces to equation 

dr dp 



r sjp 2 -2Mp 

which is simply integrated: 

ln(Cr) = In (y/ ' p 2 - 2Mp + p-M 

where C is an integration constant. Thus the transformation to the new coordinate system 
is given by the function 

Cr = vV - 2Mp + p-M. (21) 
The inverse transformation of the radial coordinate is 

Now the spatial part of the Schwarzschild metric can be written in Weyl Euclidean form: 

dl* = [dr 2 + r\dd 2 + sin W)] • 

We require this metric to be asymptotically Euclidean for r — > oo. This requirement 
uniquely determines the integration constant: C — 2. 

Finally, we obtain the following expression for the spatial part of the metric 



dV 



( x + f:) + r W + sin W)] • (23) 



We have g^ u = / 2 5 M „ in notations of the previous section. Therefore for the function / 
we get exactly expression (fTTT) . 
It is easily calculated, 

2M I 1 - 



1 - 



M N 2 
1r 



2r 

Thus the Schwarzschild metric in isotropic coordinates has the form 



ds 2 = I \ d f _ + ^ 4 [ dr * + r 2 (dd 2 + sin 2 6d V 2 )] . (24) 

For M > 0, this metric is defined everywhere in IR 4 excluding the world line of the origin 
of the spherical coordinate system: 

-oo <t<oo, 0<r<oo, 0<#<7r, < <p <2n. 

This metric satisfies vacuum Einstein's equations everywhere in M 4 excluding the world 
line of the origin of the spherical coordinate system because it is obtained by the coordinate 
transformation (|22|) from the Schwarzschild metric. According to theorem I2.1[ the space 
part of this metric is the fundamental solution of equation (jHJ) with 5-type right hand side 
for all moments of time t. 
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Metric ( 124)) is isotropic on all sections of the space-time t = const. Therefore these 
coordinates for the Schwarzschild solution are called isotropic. They are well known for 
a long time, see, i.e. [9]. The new result is not the Schwarzschild metric in isotropic 
coordinates (I2~4"l) . but the statement that its spatial part satisfies equation (j^D with S-type 
source. This fact has principal meaning and far reaching physical consequences. We shall 
see in what follows that gravitational attraction to the point mass M at large distances 
is replaced by repulsion at distances r < M/2. 

The determinant of metric (l24"j) is equal to 



det 9 „ fl = -(l-^V fl + ^Y r 4 sin 2 9. (25) 



j 2 (i+-V 

2r J V 2r J 

Consequently, the metric is degenerate on the sphere of radius 

r* := y, (26) 

which corresponds to the horizon p = 2M for the Schwarzschild solution and on the z-axis 
(6 = 0, 7r). Degeneracy on the z-axis is related to the spherical coordinate system. 

It is well known that all geometric invariants for the Schwarzschild metric constructed 
from the curvature tensor have no singularities on the horizon. 

When radius r decreases from infinity to the critical value r* := M/2, the Schwarzschild 
radial coordinate decreases from oo to the horizon p s := 2M. Afterwards the radius p 
increases from p s to oo as radius r decreases from r* to zero. We see that point massive 
particle is located at that site where the Schwarzschild radial coordinate equals infinity. 
Thus two copies of the Schwarzschild metric outside horizon 2M < p < oo are mapped 
on two distinct domains in R 3 : < r < r* and r* < r < oo. On the sphere r = r* they 
are smoothly glued together. The spatial part of the metric is not degenerate here. Note 
that area of a sphere surrounding a particle tends to infinity as it gets closer to it. It 
is related to the fact that components of the spatial part of the metric (124)) diverge for 
r — >■ 0. 

One can easily calculate the asymptotic of the zeroth component of the metric 0241) 

2M 

5-00 ~ 1 , r ^ oo. 

r 

It is the same as for the Schwarzschild metric. This is not surprising because p — Y r for 
r — > oo. It means that there is no problem with Newton's gravitational law for metric 
( |24|) because it is defined by the asymptotic at r — > oo. 



4 Point mass in General Relativity 

Let us consider topologically trivial manifold M^R 4 (space-time) with Cartesian coor- 
dinates x a , a = 0, 1, 2, 3, and metric g a /3(x) of Lorentzian signature sign g a p — (H ). 

We denote the world line of a point particle by {q a {T)} where r is a parameter along the 
world line. In General Relativity, a point particle of mass M is described by the following 
action 

S = W I dx ^\ R - M I dT \l^9a^ (27) 
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where g := detg Q/ 3, R is the scalar curvature, and q := dq/dr. 

Variations of action (|27|) with respect to metric components g a p and coordinates of a 
particle q a yield Einstein's equations of motion and equations for extremals (geodesies): 

R aP - \g afi R = ~T al3 , (28) 
q a +r^ a \ x=o fq^)g aS = 0, (29) 



=9 

where 

T a ? = j=4^5{x - q) (30) 



is the particle energy-momentum tensor, r /37 ° are Christoffel's symbols for metric g a p, 
and 

5(x - q) := 5(x l - q l )5(x 2 - q 2 )5(x 3 - q 3 ) 

is the three-dimensional 5-function on sections x° = const. We assume that for particle 
trajectory q° ^ 0. 

Usually, Einstein's equations are written with lower covariant indices. We are seeking 
solutions of Einstein's equations for metric which components may have zeroes and sin- 
gularities. Therefore equations with upper contravariant and lower covariant indices can 
be inequivalent. 

We have solved equation fl3J) for three-dimensional scalar curvature with 5-function on 
the right hand side. However it is not clear at all what relation this solution has to system 
of equations ([28]) (f29|) . 

To analyze system of equations ( I28l) - (l29l . we rewrite it in the Hamiltonian form. To 
this end, the ADM parameterization of the metric is used 

/ N 2 + N p N p N v \ . . 

9o0=[ AT „ > ( 31 ) 



Np g^ u 

where g^ v is the metric on space like sections of the space-time x° = const. In this 
parameterization, 4 functions: the lapse function N(x) and shift functions N^x), are 
introduced instead of 4 components goo and go^ of the metric. Here we use notation 
N p := g^N^ where g pp is the 3 x 3- matrix which is inverse to g^ v : 

g pfl g„ = K- 

In what follows, raising of space indexes which are denoted by Greek letters from the 
middle of the alphabet, /i, u, . . . = 1, 2, 3, is performed using the inverse three-dimensional 
metric marked with a hat. We assume that all sections x° = const are space like, and 
therefore the metric g^ v is negative definite. 

Let p pv and p a be momenta conjugate to generalized coordinates g^ and q a . The 
action for a point particle is invariant with respect to reparameterization of the world 
line. To simplify formulae, we fix the gauge r = q . Then Einstein's equations ( 128]) in the 
Hamiltonian form reduce to constraint and dynamical equations. Constraints have the 
form 

H± = - (j^jv - V 2 ) ~eR+ v^+p 2 5(x-q) = 0, (32) 
Hp = -2V u p\ - P p5{x -q) = 0, (33) 
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where p 2 := g^p^Pv Constraints (1321) and (|33|) are called dynamical and kinematical, 
respectively. Dynamical Hamiltonian Einstein's equations are 

2N N 

9fMu = —Pfiu - -jd^uP + ^^N u + V u Np, (34) 

- eN (rT - \g» v R} - P » P V P N» - p^V p N» + %(N>jr) - ^0==5(x - q), 

(35) 

where A := V^V^ is the three-dimensional Laplace-Beltrami operator and p^ := g^ v p u - 
Hamiltonian equations for geodesies have the form 



\JM 2 + p 2 



P»-N»\ x=q , (36) 



x=q 



Pn = -d P N^M 2 +p 2 -N u Pu . (37) 



x=q 



Transition from Lagrangian equations of motion ( |28i) . ( |29l) to Hamiltonian (I32 i) -(13~T I) 
is complicated. Hamiltonian formulation of General Relativity was given in [101 H]. F° r 
General Relativity and point particle it can be found i.e. in [llj . Combined Hamiltonian 
system of equations of motion for a point particle in General Relativity was considered in 

mm. 

We are seeking spherically symmetric static solution of the system of equations of 
motion ( l32l) - (l3Tj) . In this case, shift functions are zero, = 0. We choose the Weyl flat 
gauge for the spatial metric 

g^v / 

Assume also that the particle is located at the origin of the coordinate system: 

q = x° = T, <f = 0, /i = l,2,3. 

In the static case, all momenta vanish: p^ u = and p p = 0. Derivatives g pv = and 
= are also zero. Thus there are two unknown functions in this case which depend 

only on radius in spherical coordinates: N(r) and f(r). 

If the above assumptions are fulfilled then equations ( 1331) . f |34l) . and ( )36l) are identically 

satisfied. Equations ( 1321 . ( |35i) . and ( |37j) take the form 

-eR + 16irM5{x) = 0, (38) 
e (hNg^ - V^Nj - eN (h^ - ^g^Rj = 0, (39) 

d,N\ x=Q = 0, (40) 

where equation ( 137|) we divided by M. 

It is important that equation (1391) can not be divided by e, and indices /i, v can not 
be lowered because metric have singularity at r = 0, and the Ricci tensor and scalar 
curvature contain the 5-function. 
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Equation (J38]) was solved in section [2j The solution is given by equation (ITT]) . The 
difference in the sign in equation (!38|) arises because the spatial metric g^ u is now negative 
definite, and changing sign of the metric results in changing sign of the scalar curvature. 

From the Schwarzschild metric in the form (I24p . we deduce that solution of equation 
for the lapse function is 

i _ M 

N = Trb- (41) 

The lapse function is positive for r > M/2 and negative for < r < M/2. It has finite 
limit at zero, 

limiV = -1, 

r- >0 

and is infinitely differentiable for r > 0. Therefore the lapse function is locally summable 
and lies in V'(R 3 ). 

Theorem 4.1. Lapse function Hjty satisfies equation f3W in R 3 in a generalized sense. 

Proof. One can easily verify that the Laplacian of the lapse function for r > is equal to 
zero, 

AN = jd^e^duN) =0. 

Function N is continuous, and its derivatives have discontinuities at zero. Since N e 
T?'(1R 3 ) then the weak derivatives are defined everywhere, and equality 

AN = 0, 

is fulfilled in a generalized sense because the lapse function has finite limit at zero. This 
equality can be multiplied by eg^ v because near zero 

e ~ r~\ ~ r 4 . (42) 

Thus the first term on the left hand side of equation fl39l) is equal to zero in a generalized 
sense. 

One can verify that for r > the following equality holds 

eWV'N + eN (k^ - ^ 9ltl ,R) = 0. 

Hence we are left only with 5-functions which are contained in the Ricci tensor and scalar 
curvature. They do not cancel, and equation (139]) takes the form 

^-gTNeR = 4xMgTN8(x) = 0. 

This equation is fulfilled because 

g ^N ~ S^r\ r^ix) = 
near zero. □ 
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Thus we proved that the Schwarzschild metric in isotropic coordinates satisfies 
Einstein's equations (128]) for point massive particle in a generalized sense. 
We are left only with equation (HOI) for the point particle. Since 

M 1 

d r N -- 



r2 (1 + f) 2 ' 

it is not satisfied. At this point, we invoke physical arguments. The right hand side 
of equation (13"T|) is the force acting on the particle from the gravitational field which is 
produced by the particle itself. Therefore equation (137|) is understood after averaging over 
a sphere surrounding the particle. Then it is fulfilled because of the spherical symmetry. 
In other words, the particle's own gravitational field does not act on it. 

The same situation happens in Newton's gravity and electrodynamics. For example, 
for gravitational potential if = — 1/r in Newton's theory, the following equality holds 

_ 1 x" 
-0/1- = -§■ , 

and equations of motion are not satisfied. Nevertheless they are fulfilled after averaging 
over a sphere. 



5 Geodesies 

To give physical meaning for metric ( 1241) . it is necessary to analyze the behavior of 
geodesies or extremals (lines of extremal length). We assume in General Relativity that 
test particles move along geodesies. To this end we compute Christoffel's symbols 
Straightforward calculations show that only nine independent components differ from 
zero: 

r P !-g M 1 M 

°° (1 + f) 7 - 2 ' * (l + i)(l-f) ^ 



1 M 1 — — 1 

1 rr 1 _ M 2 ' 1 r0 1 r f i , M r ' 



2r 1 2r 

M 
2r 



i-^ „ , n cose 



1 + M r ' Le ^ S in0' 

I — M. 



IV = -- frsin 2 e, r w e = - sin^ cosfl. 

^ 2r 

Therefore equations for geodesies {x a (r)}, a = 0, 1, 2, 3, 



(43) 
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where x a := dx a /dr is the differentiation with respect to canonical parameter take the 
form 

2 M ■ 

1 _ M i /if i _ M 

^ _ _. 2r JKJ ^2 , 1 1VI ^2 , 1 2r 



1 ' 2r 1 ' 2r 



1 _ m 2 



(44) 



= — 2 {r -r6 + sin 6 cos #0 2 , 

2r 
M 



1-f 2.. o cos0- 
1 + t sin 6/ 

2r 

Equations for geodesies always have the integral of motion 

M \ 2 / ,/\4 



Co = x^g a , = i 2 fL-f j "(l + f) [r 2 + r 2 (£ 2 +sin 2 ^ 2 )]. (45) 

Its geometric meaning is that the length of a tangent vector is conserved along geodesic. 
Moreover, each Killing vector K = K a d a yields the integral of motion 

C = K a x l3 g a/ s = const. 

The obtained solution is static which means that K\ = d$ is the Killing vector. There 
are three Killing vectors corresponding to spherical symmetry from which only two are 
independent. Thus there are three integrals of motion for geodesies in metric (1241) related 
to the symmetry of the metric. Existence of four integrals of motion for system of equa- 
tions ( j44l) allows to analyze qualitatively the behavior of all geodesies. We analyze only 
the simplest. 

For Killing vector fields K\ = do and K 2 = d v , we have integrals of motion: 




n = *{Ttt) • (46) 

C 2 = -ip (l + ^ r 2 sin 2 #. (47) 

Differentiating Co,i,2 with respect to canonical parameter and using equations of motion 
(1441) one can prove that Co ; i >2 are integrals of motion, indeed. 

Using the freedom in the choice of a canonical parameter, we put C\ = 1 for r > r* := 
M/2. Then at large distances, r 3> M, where the gravitational field is small conservation 
law ( 14"6j) takes the form t = 1 at zeroth order. This means that at large distances the 
canonical parameter r coincides with observational time t. This gives physical meaning 
to constant C\. 

At large distances in the main approximation with respect M/r and for t — 1, second 
equation for radius in (1441) reduces to 

r = + — r 2 + r(6 2 + sin 2 6>y2 ). 
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If we require in addition smallness of velocities f < 0, then the second term on the right 
hand side can be neglected as compared with the first term. Thus at large distances and 
small velocities, 

M 



— ~ e, r ~ e, e <C 1, 
r 



equations for geodesies in General Relativity can be written in the form 

r = -^L +r (6 2 + sin 2 fly? 2 ), (4f 



r 
IrQ 



+ sin cos 6^ , (49) 



2fo9 cos 6* • , , 

(p = 1 - 2—Oip, 50 

r sin 6 

first equation in ( 1441) being satisfied in the main approximation. In this form, equations 
for geodesies can be compared with equations of motion for a point particle in classical 
Newton's theory. 

Newtonian limit. In Newton's gravity theory, a point particle of unit mass moves 
in three-dimensional Euclidean space in the potential field U = —M/r. Its motion along 
trajectory {x M (t)} is defined by the following action functional 

S= fdtL := [dt('4 + —)= fdt (rl±I^t±^lfl + —'], (51) 



2 r J J \ 2 r 

where v is the velocity of the particle which is written in spherical coordinate system. 
Variation of this functional with respect to r(t), 9(t), and <p(t) results exactly to equations 
of motion (j4"5]) - (l5"0~I) where dot denotes differentiation with respect to time t. 

Invariance of action (15T|) with respect to time translation t\-tt+ const, according to 
the Noether theorem, leads to the energy conservation 

2 M p + r 202 + sjn 2^2) M 

E = = — — — i — ILL = const. 52 

2 r 2 r y ' 

Invariance of action ( 15 ip with respect to angle shift (p i— > (p+ const leads to conservation 
of the angular momentum projection on z axis: 

L = — ^ = — r 2 sin 2 6ip = const. (53) 
dip 

Similarly, other angular momentum projections on x and y axes are conserved though 
their expressions in spherical coordinates look more complicated. 

Since the angular momentum is conserved, a particle moves on the plane which goes 
through the origin and the vector of initial velocity. Therefore without loss of generality, 
we can restrict motion to a single plane, say 9 = ir/2. 

Thus at distances much larger then r* which correspond to the Schwarzschild radius 
p s := 2M and velocities small as compared to the velocity of light, General Relativity 
reproduces equations of motion for a point particle in Newton's theory. 

Note that the Schwarzschild radius is small. Substitution of the numerical values 
yields the following values: for the Sun and Earth the Schwarzschild radii are equal to 3 
km and 1 cm, respectively. 
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Comparison of conservation laws ( 155]) and ( 152]) . and also ( H7|) and (|53l) tells us that 
integration constants Co and C2 for geodesic equations are analogs of energy and angular 
momentum conservation in General Relativity, respectively. □ 

As the consequence of angular momentum conservation (C2) we restrict ourselves to 
geodesies in the plane 9 = n/2. Then 6 = 0, and equations (jUJ) are simplified: 

2 M . 

2 tT > 



(1 + f) (1-f) 



1-f M. 2 1 M. 2 1-f . 2 

77-^ ^ + TTM ^ + T^t r ^ ■ (54) 

I 1 + IFJ 2r 1 T 2r 



1-M 2 



Integrals of motion (155j) and (147p become slightly simpler: 

^o = t 2 (^|) -(l + |) 4 [r 2 + rV], (55) 



M\ 2 



^2 = -0 ( 1 + ^ ) r". (56) 



Integral of motion (HBjl remains the same. 

For further analysis of geodesies, we put C\ = 1. This means that for infinitely distant 
observer the canonical parameter along geodesies coincides with the observed time r t. 

Circular geodesies. If a test particle moves along a circle r = const, then first 
equation in ( 15^1) takes the form t = 0. Therefore without loss of generality we can identify 
the observed time with canonical parameter, t = r. Then second equation in ( 1541) defines 
the angular velocity: 

•2 1 M 

ip = — = const. 57 

(1 + f) 6 ^ 3 

At large distances, r 3> M, the obtained angular velocity of rotations coincides with that 
in Newton's mechanics: 

ip 2 - M 

Third equation in ( |54l) for circular geodesies is automatically fulfilled. 

As the radius r decreases from infinity to the critical value r* the angular velocity of 
rotations increases from zero to its maximal value at r = r*. Then the angular velocity 
of rotations starts decreasing which means that attraction changes to repulsion. Near the 
particle it decreases to zero 

lim0 2 = 0. 

This means that near the origin of the coordinate system a test particle does not feel any 
gravitational force. 

Lightlike radial geodesies. For lightlike geodesies, Co = 0. Besides, for radial 
motion of light, = 0. In this case, conservation law ( 155]) allows to find 

(1 + K) 3 

i = ± \_V r. (58) 

2r 
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On the other hand, conservation law (HST) for C\ = 1 leads to equality 




(59) 

The obtained equations allow to find the velocity 

' = ±77^- (60) 

These equations are easily integrated: 

M 2 

± r + const = r H . (61) 

4r 

As the consequence, the points r = and r = oo are complete for radial lightlike geodesies. 
That is radially falling light never reaches a point particle of mass M. This is what is to 
be expected because the domain near the particle r <C M is diffeomorphic to the external 
Schwarzschild solution for p ^> M which is complete as p — > oo, this limit corresponding 
to r — > 0. In particular, no timelike curve can reach the origin of the coordinate system 
r = at finite proper time. 

Radial equation (154|) with equalities (1591) and fl60l yield equation for the radius 

1 M 2 

r = o— 62 

(i-IS) 32 ^ 3 

where acceleration stands on the left hand side. At large distances r > r* acceleration 
is negative, and this corresponds to attraction. At small distances r < r* acceleration is 
positive, and this means repulsion. This consequence of Einstein's equations has nothing 
in common with Newtonian gravity where we have only attraction. 

At large distances as well as near the critical radius and the origin of the coordinate 
system velocity and acceleration have the following behavior: 

r — > oo : r — >• ±1, r — > 0, 

r — > r* : r — >• ±oo, r — >■ Too, 

r — ^ : r — > 0, f — )■ 0. 

That is, at infinity and near the origin, light moves without acceleration. When crossing 
the critical sphere r = r* the acceleration changes its sign which means that attraction is 
replaced by repulsion. 



6 Conclusion 

In the present paper we prove that the Schwarzschild in isotropic coordinates satisfies 
the system of Einstein's and geodesic equations for a point massive particle. Nontrivial 
energy-momentum tensor appears on the right hand side of Einstein's equations which 
is proportional to three-dimensional 5-function. Rewriting the Schwarzschild solution in 
isotropic coordinates is by itself a student's problem. The result of the paper is that we 
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attribute mathematical meaning for Einstein's equations with 5-type source and prove 
that the Schwarzschild solution in isotropic coordinates describes gravitational field of a 
point massive particle. 

The mass M appears in the Schwarzschild solution as an integration constant and may 
be arbitrary. In our approach, the mass M of a point particle enters the action from the 
very beginning. 

The obtained solution is isometric to Einstein-Rosen bridge [13] . In the original paper 
the bridge was attributed to an elementary particle which in our notations correspond to 
the critical sphere r = r*. The space-time is described by two sheets. We have shown 
that the particle described by a (^-function is located not at r = r*, but at geodesically 
complete "infinity" of one of the sheets corresponding to r = 0. Geodesically complete 
infinity of the second sheet lies at infinity r = oo and corresponds to asymptotically 
flat space-time. Both sheets are glued together at the critical sphere r = r* where metric 
becomes degenerate and gravitational attraction at large distances is replaced by repulsion 
at small distances. This effect is of primary importance and makes us to reconsider our 
approach to many gravitational phenomena. 

At present, the Einstein-Rosen bridge has another interpretation. Two sheets are 
considered as two different universes which are connected by a worm hole at r = r*. It 
was noted recently [14] that nontrivial energy-momentum tensor appears on the right 
hand side of Einstein's equations at r#. In our approach, the right hand side of Einstein's 
equations is equal to zero at r — r*. The difference comes from choosing the solution 
to equation iV 2 = g 00 . For a given metric we chose solution (I4ip which is positive for 
r > r* and negative at r < r*. The other solution |iV| is chosen in [13], the modulus 
sign resulting in the appearance of the singularity on the right hand side of Einstein's 
equations at r = r*. 

The obtained space-time is topologically four- dimensional Euclidean space with re- 
moved world line of the particle r = 0. This space-time is geodesically complete. 

The effect of the transformation of gravitational attraction into repulsion is the straight- 
forward consequence of Einstein's equations and leads to many questions. For example, 
how can a black hole be formed if particles cannot come close one to another ? One can 
make a lot of speculations at this point but we do not consider them here. 

The author is grateful to I. V. Volovich, A. K. Gushchin, Yu. N. Drozhzhinov, B. I. Za- 
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